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The invariant differential cross section for inclusive neutral pion production in p+p collisions at 
\/s = 200 GeV has been measured at mid-rapidity (|r7| <0.35) over the range 1 <pt^ 14 GeV/cby the 
PHENIX experiment at RHIC. Predictions of next-to- leading order perturbative QCD calculations 
are consistent with these measurements. The precision of our result is sufficient to differentiate 
between prevailing gluon-to-pion fragmentation functions. 

PACS numbers: 13.85.Ni, 13.87.Fh, 25.75-q, 25.75.Dw 



Particle production at large transverse momenta, px, 
in hadronic reactions has provided an important test- 
ing ground for perturbative Quantum Chromodynam- 
ics (pQCD) [1]. Next-to-leading order pQCD calcula- 
tions describe Tevatron (-^8=1. 8 TeV) measurements of 
inclusive jet production [2] within 10% and direct pho- 



ton production [3], in which the elementary quark-gluon 
scattering produces a photon directly, within 50%. For 
high-pT pion production, the recent calculations have 
not been compared with the UA2 data [4] obtained at 
y/s=MO GeV. However, at lower center of mass energies 
(v^<63 GeV), they underestimate the data by a factor 
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of ^2.5 [5]. Similar discrepancies have been observed 
for direct photon measurements from fixed target experi- 
ments [6] and have been attributed to effects of soft-gluon 
radiations beyond NLO [7], to effects of initial intrinsic 
transverse momentum, kx [8], or to experimental prob- 
lems in the difficult direct photon measurements [9,10]. 
The TT° calculations, as compared to the jet or direct 
photon calculations, also require the probability for the 
scattered quark or gluon to fragment into a pion. 

Information on fragmentation to pions [11-15] has 
principally come from global analyses of inclusive hadron 
production in e+e~ annihilation. These analyses con- 
strain the quark-to-pion fragmentation functions well 
but, via the scale dependence, the gluon-to-pion frag- 
mentation function to a lesser extent. For example, the 
latter function at a scale of 100 GcV^ can vary by a factor 
of 2 to 10 when the fraction of the initial gluon momen- 
tum carried by the pion is above 0.5. The more direct 
measurements of gluon fragmentation functions from b- 
tagged, three jet event data from LEP [16] have played a 
limited role in the global analyses because NLO correc- 
tions are unavailable for the quantitative treatment, in- 
cluding scale and scheme dependences, of these data. As 
has been explored for measurements of inclusive hadron 
production in p+p collisions [17], results from inclusive 
pion production at high px can be included in the global 
analyses and thus may provide meaningful constraints on 
the gluon-to-pion fragmentation. These results will also 
provide a reference needed for quantifying the suppres- 
sion of 7r° production observed in Au-Au collisions at the 
Relativistic Heavy Ion Collider (RHIC) [18] and, to the 
extent of agreement with the calculations, the founda- 
tions for the planned polarized gluon density measure- 
ment with polarized protons in RHIC [19]. 

In this Letter, we report the first results on inclusive 
neutral pion production in p + p collisions at a center 
of mass energy (\/s) of 200 GeV as extracted from the 
data collected during the 2001-2002 run period (Run- 
2) of RHIC. The bunched proton beams in the collider 
were vertically polarized with spin orientations alternat- 
ing in successive bunches. By balancing the integrated 
luminosity in the different spin states, the effects from 
polarization were canceled at the 0.1% level. 

In Run-2, the PHENIX experiment [20] operated two 
central arm spectrometers, one muon arm spectrometer, 
and other detectors for triggering and vertex determina- 
tion. This work used the beam-beam counters (BBC) [21] 
for determining the collision vertex and constructing the 
minimum bias (MB) trigger, and the electromagnetic 
calorimeters (EMCal) [22] for detecting the neutral pi- 
ons and deriving high-pr triggers. 

The unbiased differential cross section for tt° produc- 
tion is calculated from the MB triggered data sample as 

^C^^f _ 1 1 C'reco • C'conv -^tt" Q^ 



where N^^o is the number of 7r°'s observed in a ApT wide 
bin at p^ defined as the pr for which the cross section 
equals its average over the bin; Ay is the rapidity range; 
Cicco is a correction for the acceptance, reconstruction 
efficiency, and pr smearing; Cconv is a correction for the 
conversion of decay photons; /tt" is the fraction of the in- 
clusive TT° yield for which the MB trigger condition was 
satisfied; and C is the integrated luminosity for the an- 
alyzed data sample. The high-pr triggered sample re- 
quired an additional correction to account for the effi- 
ciency of this trigger for 7r° detection. 

The MB trigger imposed the requirement that the col- 
lision vertex was within 75 cm of the center of the in- 
teraction region. This vertex was reconstructed from the 
difference in the arrival times of particles at the BBCs 
which were located along the beam line at ±1.44 m from 
the nominal interaction point and subtended the pseudo- 
rapidity range ±(3.0-3.9) with full azimuthal coverage. In 
the analysis of the data, a more restrictive requirement 
of ±30 cm was applied. 

The EMCal consisted of two subsystems: a six sector, 
lead scintillator (PbSc) calorimeter and a two sector, lead 
glass (PbGl) calorimeter. Located at a radial distance of 
~5 m from the beam line, each of these sectors covered 
the pseudorapidity range of |r?|<0.35 and an azimuthal 
angle interval of A(/)w22.5°. Each of the towers in the 
calorimeter subtended A0x A?^ ~ 0.01x0.01, thus ensur- 
ing that the two photons from a decayed 7r° were resolved 
up to a pt of at least 20 GcV/c. The energy calibration 
was corroborated by the position of the 7r° invariant mass 
peak, the energy deposit from minimum ionizing charged 
particles traversing the EMCal (PbSc), and the correla- 
tion between the energy deposit in the EMCal and the 
measured momentum for electrons and positrons identi- 
fied by the ring-imaging Ccrcnkov detector. These stiid- 
ies showed that the accuracy of the energy measurement 
was within 1.5%. At a pr of -^11 GeV/c, this uncer- 
tainty translates into a systematic error on the tt° yield 
of ~12%. The effective energy resolution for the dataset 
was deduced from the widths of the 7r° mass peaks, which 
varied with PT from 7% to 10% (PbSc) and 12% to 13% 
(PbGl), and a comparison of the measured energy and 
momentum for the identified electrons and positrons. 

The number of recorded high-pT 7r°'s was enhanced 
by a high-pT trigger (denoted as 2x2) in which thresh- 
old discrimination was applied independently to sums 
of the analog signals from non-overlapping, 2x2 group- 
ings (called tiles) of adjacent EMCal towers. During this 
run, the thresholds corresponded to a deposited energy 
of 0.75 GeV. The efficiency of this trigger for 7r° de- 
tection, £^o^(pt), was obtained from the MB data. As 
shown in Fig. la, this efficiency reached a plateau at a 
Pt of ~3 GeV/c. This despondence was reproduced by 
Monte Carlo calculations which included the measured 
tile threshold curves, the EMCal detector response, and 



3 



the geometry of the active trigger tiles. The saturation 
level, 0.78 ± 0.03 for both PbSc and PbGl, was consis- 
tent with the geometrical acceptance of the active trigger 
tiles. For pT>4: GeV/c, the geometrical acceptance was 
used to correct the tt° yield in the high-pT sample. 

Since only a fraction of inelastic p+p collisions produce 
particles which enter both BBCs, the MB trigger condi- 
tion biased the recorded data sample, so only a fraction, 
/^o, of the inclusive tt° yield was detected. This frac- 
tion was determined from data collected by an additional, 
high-pT trigger which had not been operated in coinci- 
dence with the MB trigger. This trigger was formed by 
threshold discrimination of the sums of the analog signals 
from overlapping 2x2 groupings of adjacent 2x2 trigger 
tiles. As shown in Fig. lb, the fraction of these high-pr 
events with 7r°'s for which the MB condition was also 
satisfied was 0.75 ± 0.02, independent of pr- 

In each event, the two photon invariant mass was cal- 
culated for each pairing of clusters. Clusters were paired 
if the energy asymmetry, \Ei — E2\/{Ei + E2), was less 
than 0.8 (PbSc) or 0.7 (PbGl). For the PbGl, the pair- 
ings were further restricted to those clusters which were 
identified as electromagnetic via the shower profile and 
time-of-flight. The 7r° yield was extracted by integrating 
the invariant mass spectrum over a region around the 
7r° mass. The background contribution in each px bin 
was estimated and then subtracted by fitting the invari- 
ant mass distribution outside the peak region (PbSc) or 
using the mixed event technique (PbGl). For the PbSc 
and the PbGl, the background to signal ratio varied with 
increasing pT from 1 to 0.1 and 1 to 0.03, respectively. 

The raw yields were corrected for the pt smearing aris- 
ing from the EMCal energy resolution and the steeply 
falling spectrum; and for the losses due to the disabled 
towers, the incomplete azimuthal coverage, the energy 
asymmetry cut, and the photon identification cut (PbGl). 
The correction for these effects, Cicco, was calculated 
with Monte Carlo siirmlations which contained the con- 
figuration of the active EMCal towers. The energy and 
position of the decay photons were smeared with the mea- 
sured test beam resolutions [22] augmented by a constant 
energy smearing of 5% (12%) for the PbSc (PbGl) to 
match the response of the EMCal. 

The correction for the losses due to conversions of de- 
cay photons, Cconv, was determined to be 4% (PbSc) 
and 9% (PbGl) by using a GEANT3 [23] simulation of 
the PHENIX detector. The same simulation, using p+p 
events from the PYTHIA generator [24], showed that 
the contribution of tt° 's from secondary interactions was 
negligible and that the contribution from decays of other 
hadrons {e.g., K° and r] mesons) was less than 6%. The 
7r° spectrum was not corrected for these decays. 

The integrated luminosity, £, was determined from the 
number of MB events using an absolute calibration of the 
trigger cross section obtained via the van der Meer scan 
technique [25]. In a scan, the transverse profile of the 



beam overlap is measured by sweeping one beam across 
the other in steps while monitoring the MB trigger rate. 
This information, the bunch intensities of the two beams 
(^lO'^-'^/bunch), and the revolution frequency (78 kHz) 
are then used to compute the instantaneous luminosity. 
The trigger cross section is the ratio of the MB trig- 
ger rate when the beams were overlapping maximally to 
the instantaneous luminosity. Based on three scans, this 
cross section was 21.8 ±0.9 (2.8) mb at the 68.5% (95%) 
confidence level with an absolute error of 0.7 mb. From 
the linear sum of the absolute error and half of the 95% 
confidence level, point-to-point systematic error, an error 
of 9.6% was assigned for the luminosity normalization. 

During the run, the maximum and average instanta- 
neous luminosities were 1.5x10^° and 0.5x10^° cm~^-s~^, 
respectively. Contributions from multiple collisions per 
bunch crossing and beam-gas interactions were negligible. 
The MB trigger sample of 16 million events corresponded 
to 0.7 nb~^. As computed from the fraction of recorded 
MB events which also met the 2x2 high-pT trigger con- 
dition (^1/47), the 18 million high-pr triggered events 
corresponded to an effective luminosity of 39 nb"""^. 

The invariant differential cross sections obtained from 
the MB and high-pT samples were consistent within 
the statistical errors over the px region of overlap 
(pt<5.5 GeV/c). Moreover, the results determined in- 
dependently from the PbSc and the PbGl data samples 
were consistent. The main sources of the point-to-point 
systematic uncertainty in the two measurements are sum- 
marized in Table I for a low and a high pt bin. The total 
error was computed as the quadrature sum of the statis- 
tical and point-to-point systematic errors. 

From the MB and the high-p^ trigger samples for px 
below and above 4 GeV/c, respectively. Table 11 tabu- 
lates the cross section and the errors obtained by aver- 
aging the PbSc and PbGl results using the total error 
for the weighting. Figs. 2a and 2b show this combined 
result and its fractional statistical and systematic uncer- 
tainties (Act/ct). The data are well parameterized by 
a power-law form A - (l+p^/po)^" with parameters of 
yl=386 mb-GeV-2.c'\ po=1.219 GeV/c, and n=9.99. 

In Fig. 2, our results are compared with NLO pQCD 
calculations [26-28] . The basic concept underlying these 
calculations is the factorization of the cross section into 
parton distribution functions for the protons, parton- 
to-pion fragmentation functions, and short-distance par- 
tonic hard-scattering cross sections which can be evalu- 
ated using perturbative QCD. Because of this factoriza- 
tion, the calculations depend on unphysical, factorization 
and renormalization scales which are of the order of the 
hard scale pr- This dependence is reduced as higher or- 
der terms are included in the perturbation expansion. 
For a calculation truncated at a given order, this depen- 
dence serves as a gauge for the uncertainty in its results. 

The calculations in Fig. 2 have been performed with 
equal renormalization and factorization scales of Pt/2, 
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PT, and 2pT by using the CTEQ6M [29] set of parton 
distribution functions and two sets of fragmentation func- 
tions. In general, these calculations are consistent with 
the data, even at low pr (<2 GcV/c) where the theory 
might be expected to be less applicable. On closer in- 
spection, as shown Fig. 2c and 2d, the calculation with 
the "Kniehl-Kramcr-Potter" (KKP) set of fragmentation 
functions [12] agrees with our data better than the cal- 
culation with the "Kretzer" set [13] does. These two 
sets differ mainly in that the gluon-to-pion fragmentation 
function, , is greater in the KKP set. This difference is 
exhibited primarily at low pt because of the dominance 
of gluon-gluon and gluon-quark interactions for pT below 
~10 GeV/c [26]. Our measurement thus may impose a 
meaningful constraint on . 

In summary, the invariant differential cross section for 
inclusive neutral pion production m p-^p collisions at 
^/s = 2Q{) GeV was measured at mid-rapidity (jr/j <0.35) 
as a function of pr up to ~14 GeV/c. These results were 
compared with two NLO pQCD calculations which dif- 
fered in the choice of fragmentation functions. Over the 
full range of py, the calculations were consistent with 
the result within the uncertainty of the calculations as 
judged from the scale dependence, although the results 
favored a larger gluon-to-pion fragmentation function. 

We thank the staff of the Collider- Accelerator and 
Physics Departments at BNL for their vital contribu- 
tions. We thank Werner Vogelsang and Stefan Kret- 
zer for their interest and input. We acknowledge sup- 
port from the Department of Energy and NSF (U.S.A.), 
MEXT and JSPS (Japan), CNPq and FAPESP (Brazil), 
NSFC (China), CNRS-IN2P3 and CEA (France), BMBF, 
DAAD, and AvH (Germany), OTKA (Hungary), DAE 
and DST (India), ISF (Israel), KRF and CHEP (Ko- 
rea), RAS, RMAE, and RMS (Russia), VR and KAW 
(Sweden), U.S. CRDF for the FSU, US-Hungarian NSF- 
OTKA-MTA, and US-Israel BSF. 



* Deceased 

^ Spokesperson: zajc@nevis.columbia.edu 
[1] J. F. Owens, Rev. Mod. Phys. 59 465 (1987). 
[2] D. Stump et al, hep-ph/0303013. 

[3] D. Acosta et al. [CDF Collaboration], Phys. Rev. D65, 
112003 (2002); V.M. Abazov et al. [DO Collaboration], 
Phys. Rev. Lett. 87, 251805 (2001). 

[4] M. Banner et al. [UA2 Collaboration], Phys. Lett. B115 
59 (1982). 

[5] P. Aurenche et al, Eur. Phys. J. C13, 347 (2000). 

[6] W. Vogelsang and M.R. Whalley, J. Phys. G: Nucl. Part. 

Phys. 23, Al (1997). 
[7] E. Laenen et al, Phys. Rev. D63, 114018 (2001). 
[8] L. Apanasevich et al, hep-ex/0204031. 



[9] P. Aurenche et al, Eur. Phys. J. C9, 107 (1999). 
[10] L. Apanasevich et al, Phys. Rev. D63, 014009 (2001). 
[11] J. Binncwics, B.A. Kniehl, and G. Kramer, Phys. 

Rev. D52, 4947 (1995). 
[12] B.A. Kniehl et al, Nucl. Phys. B597, 337 (2001). 
[13] S. Krctzcr, Phys. Rev. D62, 054001 (2000). 
[14] S. Kretzer et al, Eur. Phys. J. C22, 269 (2001). 
[15] O. Biebel et al, hep-ph/0109282. 

[16] P. Abreu et al [DELPHI Collaboration], Eur. Phys. 
J. C17, 207 (2000); G. Abbiendi et al [OPAL Collab- 
oration], Eur. Phys. J. Cll, 217 (1999); R. Baratc et al 
[ALEPH Collaboration], Eur. Phys. J. C17, 1 (2000). 

[17] L. Bourhis et al, Eur. Phys. J. C19 89 (2001). 

[18] S.S. Adler et al [PHENIX Collaboration], submitted to 
Phys. Rev. Lett. 

[19] G. Bunce et al, Ann. Rev. Nucl. Part. Sci. 50, 525 
(2000). 

[20] K. Adcox et al, Nucl. Instrum. Methods in Phys. Re- 
search A499, 469 (2003). 

[21] M. Allen et al, Nucl. Instrum. Methods in Phys. Re- 
search A499, 549 (2003). 

[22] L. Aphecetche et al, Nucl. Instrum. Methods in Phys. 
Research A499, 521 (2003). 

[23] GEANT 3.21, CERN program library 

[24] T. Sjostrand, Comp. Phys. Commun. 82, 74 (1994). 

[25] K. A. Drees and Z. Xu, "Proceedings of the PAC2001 
Conference," 3120 (2001). 

[26] F. Aversa et al, Nucl. Phys. B327, 105 (1989); n.h., 
these authors wrote the computer code used for our cal- 
culations. 

[27] B. Jagcr et al, Phys. Rev. D67, 054005 (2003). 
[28] D. de Florian, Phys. Rev. D67, 054004 (2003). 
[29] J. Pumplin et al [CTEQ Collaboration], J. High Energy 
Phys. 0207, 012 (2002). 

TABLE I. Summary of the sources of systematic errors on 
the 7r° yields and the total systematic error for pr of 1.2 and 
10.9 GeV/c. The normalization error of 9.6% is not listed. 



% Error (PbSc) % Error (PbGl) 



Pt (in GoV/c) 


1.2 


10.9 


1.2 


10.9 


Energy Scale 


3 


11 


6 


12 


Yield Extraction 


7 


4 


5 


5 


Yield Correction 


3 


6 


6 


11 


Acceptance Stability 


4.5 


4.5 


3 


2 


Total 


9 


14 


10 


17 
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TABLE II. The (see text for definition), the invariant 
difi'erential cross section for inclusive neutral pion production 
in p+p collisions at -^s = 200 GeV, the statistical uncertainty, 
and the systematic uncertainty for each pr bin. The absolute 
normalization error of 9.6% is not included. 







inv. 


cross 


Stat. 


syst. 


Pt bin 


Pt 


section 


error 


error 


(GeV/c) 


(GeV/c) 


(mb-GeV-^-c^) 


(%) 


(%) 


1.0-1.5 


1.22 


3.73 


10-^ 


1.6 


7.3 


1.5-2.0 


1.72 


6.05 


10"^ 


1.8 


7.1 


2.0-2.5 


2.22 


1.22 


10"^ 


2.5 


7.1 


2.5-3.0 


2.73 


3.31 


10"^ 


3.6 


7.2 


3.0-3.5 


3.23 


9.98 


10-* 


5.7 


7.3 


3.5-4.0 


3.73 


3.39 


10-* 


7.3 


7.7 


4.0-4.5 


4.23 


1.19 


10-* 


2.4 


8.3 


4.5-5.0 


4.73 


4.73 


10"^ 


4.2 


8.5 


5.0-5.5 


5.23 


2.21 


10-^ 


5.0 


8.7 


5.5-6.0 


5.74 


1.11 


10"^ 


4.5 


9.2 


6.0-6.5 


6.24 


5.00 


10-^ 


6.3 


9.5 


6.5-7.0 


6.74 


3.00 


10-® 


7.7 


9.8 


7.0-8.0 


7.45 


1.08 


10-® 


8.8 


10.1 


8.0-9.0 


8.46 


4.85 


10-^ 


12.0 


10.8 


9.0-10.0 


9.46 


1.64 


10-^ 


19.3 


11.0 


10.0-12.0 


10.86 


5.07 


10-® 


22.3 


11.7 


12.0-15.0 


13.25 


9.76 


10-*^ 


41.3 


15.6 
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FIG. 1. a) The efficiency of the 2x2 high-pT trigger for 
7r° 's as a function of the pt of the 7r° . The dashed and solid 
lines show the results of a Monte Carlo simulation based on 
the 2x2 trigger tile efficiencies and the limit derived from the 
fraction of active trigger tiles, respectively, b) The fraction 
of the inclusive 7r° yield which satisfied the MB trigger con- 
dition. The solid line shows a fit of these data to a constant. 
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FIG. 2. a) The invariant differential cross section for in- 
clusive n° production (points) and the results from NLO 
pQCD calculations with equal renormalization and factor- 
ization scales of pt using the "Kniehl-Kramer-Potter" (solid 
line) and "Kretzer" (dashed line) sets of fragmentation func- 
tions, b) The relative statistical (points) and point-to-point 
systematic (band) errors. c,d) The relative difference between 
the data and the theory using KKP (c) and Kretzer (d) frag- 
mentation functions with scales of pt/2 (lower curve), pr, 
and 2pT (upper curve). In all figures, the normalization error 
of 9.6% is not shown. 
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